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DRAG  MEASUREMENTS  ON  LONG,  THIN  CYLINDERS  AT 
SMALL  ANGLES  AND  HIGH  REYNOLDS  NUMBERS 


1.  INTRODUCTION 


Historically,  the  problem  of  axisymmetric  turbulent  boundary  layers  that  develop  on 
long,  thin  cylinders  at  small  angles  of  attack  has  received  much  less  attention  than  has  the 
analogous  flat-plate  problem.  The  “axisymmetric  problem”  relates  directly  to  towed  sonar  arrays 
in  turns,  in  maneuvers  and  transient  motions,  and  in  cases  where  arrays  experience  positive  or 
negative  buoyancy  in  steady-state  tows.  Existing  modeling  and  resulting  real-time  towed  array 
shape  predictions  are  limited  in  accuracy  by  the  uncertainty  in  the  orthogonal  and  tangential  drag 
coefficient  values.  The  turbulent  wall  pressure  fluctuations,  as  well  as  the  flow-induced  cylinder 
vibrations  for  this  class  of  flows,  are  also  of  primary  importance  with  regard  to  self-noise  of 
towed  array  sonar  systems.  Willmarth  et  al.  and  Bull  and  Dekkers2  suggested  that  significant 
effects  on  the  turbulent  boundary  layer  structure  may  occur  because  of  very  small  angles  of 
attack.  This  sensitivity  to  small  angles  led  Willmarth  to  design  a  vertical  wind  tunnel  to 
eliminate  cylinder  static  deflections  caused  by  gravity. 

For  very  long  cylinders,  the  thickness  of  the  boundary  layer  (<5)  in  relation  to  the  cylinder 
radius  ( a ),  that  is,  51a,  becomes  an  additional  parameter  of  importance.  Measurements  by 
Cipolla  and  Keith3  conducted  at  the  National  Aeronautics  and  Space  Administration  (NASA) 
Langley  High-Speed  Tow  Tank,  for  which  5/ a  >  100  and  LI  a  >  10,000  (where  L  is  cylinder 
length),  established  that,  compared  to  flat-plate  cases,  axisymmetric  boundary  layers  have  a 
significantly  higher  wall  shear  stress  and  a  slower  spatial  growth  rate. 

Cipolla  and  Keith4  investigated  zero-pressure  gradient  turbulent  boundary  layers 
developed  on  cylinders  with  diameters  of  0.61,  0.89,  and  2.5  mm  and  Reynolds  numbers  of 
104  <  Re&<  105  and  108  <  ReL  <  109.  The  results  showed  that  the  wall  shear  stress  is  significantly 
higher  than  that  for  a  comparable  flat-plate  case  and  the  spatial  growth  of  the  boundary  layers  is 
lower  than  that  for  a  comparable  flat-plate  case.  Also,  the  mean  wall  shear  stress  exhibited 
spatial  variations  not  seen  in  zero-pressure  gradient  flat-plate  turbulent  boundary  layers.  For 
example,  measurements  on  towed  lines  with  radius  a  =  0.445  mm  and  L  =  137  m  have  shown 
that  the  mean  wall  shear  stress  increases  to  a  maximum  at  approximately  50  m  from  the  tow 
point,  decreases  to  a  minimum  at  90  m,  and  then  begins  to  increase  again.  These  length  scales 
are  three  orders  of  magnitude  larger  than  the  maximum  momentum  thickness  at  the  end  of  the 
line.  Assuming  the  boundary  layer  thickness  is  an  order  of  magnitude  greater  than  the 
momentum  thickness  (which  is  correct  for  flat-plate  cases),  a  length  scale  of  90  m  is  two  orders 
of  magnitude  greater  than  5  for  these  cases.  Typically,  10<5  is  considered  the  maximum  length 
over  which  boundary  layer  turbulence  displays  any  significant  correlation.  These  results  suggest 
that  a  length  scale  on  the  order  of  100 5  exists  to  observe  this  effect.  For  this  diameter 
(0.89  mm),  the  effect  was  found  for  towing  speeds  of  3.1,  9.3,  and  14.4  m/s,  with  the  minimum 
shear  stress  occurring  at  approximately  the  same  streamwise  distance  (90  m)  for  each  tow  speed. 
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Furthermore,  ongoing  stereo  particle  image  velocimetry  measurements  by  Furey  et  al.5 
for  towed  lines  of  the  same  diameter  reveal  a  decrease  in  the  boundary  layer  thickness  to  a 
minimum  value  at  the  same  location  (90  m)  as  the  minimum  in  the  shear  stress.  The  apparent 
quasi-periodic  behavior  in  both  8  and  r„  (axial  wall  shear  stress)  has  never  been  experimentally 
observed  before.  In  each  of  these  investigations,  the  tow  angle  was  less  than  1°.3'5 

In  this  investigation,  the  problem  of  a  long,  thin  cylinder  at  very  small  angles  to  the  flow 
at  moderate-to-high  Reynolds  numbers  is  considered.  This  geometry  has  several  applications  in 
which  the  total  drag  in  the  direction  of  tow  is  a  primary  measured  quantity.  In  this  research,  the 
analysis  is  extended  to  calculate  the  force  along  the  axis  of  the  cylinder,  which  is  a  primary 
design  parameter  for  towed  sonar  arrays.  This  approach  allows  an  evaluation  of  the  growing 
turbulent  boundary  layer,  which  is  particularly  relevant  for  cases  of  small  angles  of  attack.  The 
force  orthogonal  to  the  cylinder  axis  is  also  calculated  and  is  used  to  determine  a  scaling  that 
accounts  for  the  effect  of  the  turbulent  boundary  layer  in  the  axial  direction  on  this  force. 

The  use  of  a  high-speed  towing  tank  generates  length  Reynolds  numbers  and  ratios  of  the 
outer  to  inner  boundary  layer  length  scales  that  are  directly  relevant  to  undersea  towed  array 
applications.  This  investigation  required  LI  a  to  be  on  the  order  of  100  and  required  the  length  to 
be  large  in  comparison  to  relevant  turbulent  length  scales  generated.  Despite  the  large  size  of  the 
facility,  the  tank  depth  limits  the  length  of  cylinder  that  can  be  towed  at  an  angle.  Velocity 
profile  measurements  were  not  available  at  the  time  of  this  investigation.  A  direct  measurement 
of  the  total  drag  on  long,  small-diameter  cylinders  enables  very  accurate  measurements  of  the 
spatially  averaged  mean  wall  shear  stress.  This  measurement,  combined  with  measurements  of 
the  cylinder  angle,  leads  to  a  relationship  among  the  measured  drag,  axial  shear  force,  orthogonal 
force,  cylinder  angle,  and  tow  speed.  A  control  volume  analysis  was  used  to  determine  the 
momentum  thickness  6  at  the  end  of  the  cylinder  from  the  axial  force. 
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2.  DESCRIPTION  OF  EXPERIMENTS 


2.1  TEST  ENVIRONMENT 

Towing  tests  were  conducted  using  carriage  number  5  at  the  Naval  Surface  Warfare 
Center  (NSWC),  Carderock  Division,  High-Speed  Towing  Basin  (see  Schoenherr  and  Brownell6 
or  go  to  http://www50.dt.navy.mil/facilities/data/carr5.html  for  the  specifications  of  this  facility). 
The  High-Speed  Towing  Basin  comprises  two  adjoining  sections:  (1)  a  deep-water  section  that 
is  4.9  m  deep,  514  m  long,  and  6.4  m  wide  and  (2)  a  shallow- water  section  that  is  3  m  deep, 

356  m  long,  and  6.4  m  wide. 

Figure  1  is  a  cross-sectional  schematic  of  the  shallow  end  of  the  tow  basin.  The  carriage 
is  propelled  by  16  weight-bearing  vertical  drive  wheels  and  16  horizontal  driving  guide  wheels. 
The  speed  is  regulated  with  an  adjustable,  direct-current-voltage,  automatic  feedback, 
computerized  control  system.  An  automatic  controller  maintains  the  carriage  speed  to  within 
±0.03  m/s.  The  maximum  carriage  speed  is  25.7  m/s  (50  knots),  and  the  maximum  average 
acceleration  rate  is  0.49  m/s2.  Additionally,  wave-absorber  troughs,  which  dampen  standing 
surface  waves  generated  by  the  tow  strut,  extend  along  the  side  walls  of  the  tow  tank.  The 
freshwater  in  the  tank  was  at  a  uniform  temperature  of  20.6°C,  with  a  kinematic  viscosity  of  1.01 
x  10"6  m2/s  and  a  density  of  998  kg/m3.  No  measurable  change  in  the  water  density  and  salinity 
was  measured  throughout  the  depth  and  length  of  the  tank  over  the  duration  of  the  experiments. 


Figure  1.  Cross-Sectional  Schematic  of  Carriage  Number  5  and  the  Towing  Basin 
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2.2  TEST  SETUP 


The  test  setup  was  configured  for  towing  cylinders  ( L  =  3.05  and  6.10  m,  a  =  6.35  mm)  at 
critical  angles  over  a  range  of  speeds  in  a  straight  tow  configuration,  with  a  simple  pinned 
connection  at  the  tow  point.  For  a  particular  tow  speed,  variations  in  angle  were  achieved  by 
interchanging  cylinders  of  different  densities.  This  setup,  shown  in  figure  2,  eliminated  the  need 
for  support  at  the  trailing  edge,  which  would  cause  additional  unknown  drag  forces  applied  to  the 
load  cell.  This  design  also  leads  to  a  zero-bending  moment  at  the  tow  point  and  load  cell  and 
minimizes  bending  along  the  length  of  the  cylinder.  A  cutaway  view  of  the  tow  point  inside  the 
fairing  is  shown  in  figure  3. 


Figure  2.  Schematic  of  the  Towing  Configuration 


Figure  3.  Solid  Model  of  the  Tow  Point  Design 
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The  load  cell  was  attached  to  the  base  of  the  towing  strut,  and  a  linear  bearing  was  used 
to  support  a  connecting  rod  between  the  load  cell  and  the  connecting  pin.  The  fixture  was 
designed  to  minimize  transmission  of  any  vertical  force  to  the  load  cell  and  was  enclosed  in  a 
fiberglass  fairing  to  minimize  the  generation  of  turbulence.  For  each  time  series  of  load  cell 
data,  the  mean  value  was  calculated  while  the  carriage  speed  was  at  steady  state.  This 
ensemble-averaged  force  represents  the  spatially  integrated  component  of  force  in  the  tow 
direction  over  the  surface  of  the  cylinder.  For  small  angles,  this  force  is  dominated  by  the 
streamwise  component  of  the  mean  wall  shear  stress,  with  an  additional  component  caused  by 
the  pressure  and  circumferential  shear  forces  on  the  cylinder  surface.  Because  of  the  large  aspect 
ratios,  LI  a  =  480  and  960,  the  effects  of  cylinder  base  drag  were  negligible.  This  result  was 
verified  experimentally  by  repeating  drag  measurements  with  and  without  two  different  tail 
cones  attached  to  the  cylinder  trailing  edge.  The  tail  cones  were  right  circular  cones  with 
diameters  of  12.7  mm  and  lengths  of  76.2  mm  and  38.1  mm,  respectively.  The  presence  of  either 
cone  caused  no  increase  in  measured  drag  with  respect  to  the  blunt  end  case.  For  each 
combination  of  tow  speed  and  tow  angle,  multiple  runs  were  performed  to  demonstrate 
repeatability. 

Cylinder  tow  angles  were  measured  using  a  Dantec  HiSense  digital  camera,  with  an  8-bit 
resolution  and  a  1280  x  1024-pixel,  60-mm  lens.  The  camera  was  located  in  a  photography  pit, 
488  m  down  the  tow  tank.  A  photograph  of  the  aluminum  cylinder  towed  at  13  m/s  is  shown  in 
figure  4,  where  the  area  of  each  square  in  the  black  background  grid  is  0.3  m2. 


Figure  4.  Digital  Photo  of  an  Aluminum  Cylinder  Towed  at  13  m/s 


It  was  determined  from  the  digital  photographs  that  all  cylinders  were  towed  in  a  linear 
configuration,  with  a  standard  error  of  less  than  5%.  The  matrix  of  completed  test  runs  is  given 
in  table  1,  which  lists  the  cylinder  material  used  to  achieve  each  combination  of  tow  angle  and 
speed. 
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Table  1.  Tow  Angles  Versus  Tow  Speeds  for  the  Materials  Tested 


Speed  (knots) 

Measured  Angle  (degrees) 

1.5 

2 

2.5 

3 

3.5 

O 

5 

5.5 

6 

6.5 

7 

m 

11.5 

5 

10 

PVC 

AL 

12 

AL 

iH 

■ 

SS 

15 

PVC 

AL 

Ti 

SS 

20 

AL 

Ti 

25 

AL 

Ti 

SS 

30 

35 

40 

Ti 

□ 

SS  =  Stainless  Steel.  AL  =  Aluminum.  Ti  =  Titanium.  PVC  =  Polyvinyl  Chloride. 


2.3  TEST  METHODOLOGY 
2.3.1  Fluid  Forces 

In  addition  to  the  two  measured  quantities  Fx  and  a,  the  fluid  forces  acting  on  the 
cylinder  caused  by  the  turbulent  flowfield  are  depicted  in  figure  5.  As  shown,  these  forces  are 
considered  to  be  steady-state  values  averaged  over  the  length  of  the  cylinder;  the  tow  direction  is 
left  to  right.  The  resultant  forces  are  placed  at  the  cylinder’s  center  of  mass,  which  assumes  that 
the  forces  vary  linearly  along  the  length.  The  error  in  this  assumption  is  considered  negligible 
because  of  the  slow  growth  of  cylindrical  boundary  layers  found  in  previous  experiments.  For 
the  range  of  angles  considered  in  these  experiments  ( a  <  12°),  the  boundary  layer  was  assumed 
to  be  attached — but  not  symmetrical.  The  xy-coordinate  system  is  moving  at  a  constant  velocity 
equal  to  the  tow  speed.  Because  the  cylinder  is  towed  in  a  critical  angle  configuration  with  a 
pinned  connection  at  B,  MB  -  0. 


Figure  5.  Free-Body  Diagram  of  a  Cylinder  Towed  in  a  Critical  Angle  Configuration 


The  objective  of  the  following  analysis  was  to  develop  an  expression  for  the  fundamental 
boundary  layer  parameters  r„  and  <9  in  terms  of  the  measured  drag  and  angle.  Figure  6  is  a 
cross-sectional  view  of  the  cylinder,  normal  to  the  cylinder  axis,  showing  the  pressure  and 
circumferential  shear  stress  in  that  plane. 


Figure  6.  Cross-Sectional  View  of  Cylinder 


Considering  the  forces  in  the  fiy- plane,  the  components  of  pressure  force  are 

2  K  L 

Pp  =  J  ^ p(<fi)a  cos  (1) 

0  0 

and 

2k  L 

Py  =  J  ^p{<f>)a  sin  (jritfxli;.  (2) 

o  o 


The  circumferential  wall  shear  stress  also  yields  force  components  in  the  /?-  and  y- directions: 

2  K  L 

Tp  =  J  ^a{(f)a  cos  (fdtpd^,  (3) 

o  o 

and 

2k  L 

7^=1  ^(7{<f)a  sin  (fd<f>di;.  (4) 

o  o 
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Because  of  the  symmetry  of  the  flow  (zero  yaw  angle),  Pp=  0  and  Tp=  0.  The  forces 
P y  and  Ty  contribute  to  the  x-  and  y-components  of  force  shown  in  figure  5.  The  force  attributed 
to  the  wall  shear  stress  in  the  axial  direction  is 

5=  (5) 

0  0 

The  forces  given  in  equations  (2),  (4),  and  (5)  comprise  the  total  horizontal  and  vertical  fluid 
forces  acting  on  the  cylinder: 

Fx  =  Py  sina  +  rr  sin  a  +  5  cos  a,  (6) 

and 

Fy  =  Pycosa  +  Trcosa-Ssina.  (7) 

Static  equilibrium  defines  the  reaction  forces  at  the  tow  point  as 

Fx=Rx>  (8) 

and 


The  following  expressions  for  the  axial  wall  shear  stress  and  tangential  drag  coefficient  follow 
immediately  from  equation  (12): 
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(13) 


T 


W 


S 


and 


Q  = 


(14) 


This  formulation,  which  gives  the  shear  force  tangential  to  the  cylinder,  is  most  relevant 
to  the  case  of  a  cylindrical  turbulent  boundary  layer  and  enables  the  calculation  of  the 
momentum  thickness  from  a  control  volume  analysis,  which  is  analogous  to  what  is  stated  by 
Cipolla  and  Keith3  and  is  described  in  section  2.3.2.  In  addition,  the  force  along  the  axis  of  the 
cylinder  is  of  primary  interest  to  towed  array  applications,  where  the  tension  in  the  array  is  a 
critical  design  parameter.  This  force  is  also  required  for  dynamic  models  of  towed  array  shape, 
in  which  bending  forces  within  the  array  are  often  neglected.  Because  towed  arrays  are  typically 
towed  at  small  angles  of  attack,  these  experiments  are  representative  of  the  type  of  flow 
experienced  in  Navy  applications,  although  the  cylinders  described  here  have  a  much  smaller 
aspect  ratio  {LI a). 

The  total  force  exerted  by  the  fluid  on  the  cylinder  orthogonal  to  the  cylinder  axis  is 
composed  of  the  pressure  and  the  circumferential  shear  forces  in  the  /-direction: 


F  =P +T . 

Y  Y  Y 


(15) 


Combining  equations  (6)  and  (11)  leads  to  an  expression  for  {Py  +  Tr)\ 

Fy=Pr  +  Tr=  Wnet  cos  a.  (16) 

An  orthogonal  drag  coefficient  can  then  be  defined  as 

C,=-, - ^ - •  (17) 

-pin,  smaf  A, 


Note  that  this  coefficient  uses  the  relevant  component  of  velocity,  U0  sin  a,  to  normalize  the 
force. 


2.3.2  Momentum  Thickness 

In  figure  7,  a  cylindrical  control  volume  is  defined  so  that  that  the  outer  edge  is  the 
maximum  boundary  layer  thickness  at  the  end  of  the  cylinder.  The  axis  of  the  control  volume  is 
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coincident  with  the  axis  of  the  cylinder,  that  is,  the  £-axis.  Note  that  this  layer  is  an  asymmetric 
boundary  layer,  so  the  resulting  values  of  8  and  momentum  thickness  6  will  be  average  values 
over  the  circumference.  These  values  are  relevant  to  towed  arrays  because  the  internally 
mounted  hydrophones  are  not  directional  and  therefore  average  around  the  circumference. 
Furthermore,  flow  noise  models  employ  average  values  of  the  boundary  layer  parameters  at  a 
given  axial  location.  The  average  momentum  thickness  can  serve  as  an  outer-length  scale,  and 
r„.  is  required  for  the  inner-length  scale — which  are  the  relevant  length  scales  for 
boundary-layer-dominated  flows.  For  cases  of  larger  tow  angles,  for  which  the  boundary  layer 
separates,  the  cylinder  diameter  then  becomes  a  relevant  length  scale. 


Figure  7.  Control  Volume  for  a  Cylinder  at  an  Angle  to  the  Tow  Direction 


When  considering  the  forces  on  the  control  volume  in  the  ^-direction,  the  only  surface 
force  is  the  shear  force  S  acting  on  the  wall  of  the  cylinder.  Control  surfaces  1  and  2  (figure  7) 
experience  different  static  pressure  forces  caused  by  the  differences  in  elevation.  This  pressure 
differential  is  balanced  by  the  component  of  the  net  weight  in  that  direction,  so  that  the  net  body 
force  Fb  is  zero.  For  the  case  of  steady-state  flow,  the  governing  equations  (conservation  of 
mass  and  momentum)  are 

jpF*dA=0,  (18) 

cs 


and 


\u(pV  •dA=S  +  FB{.  (19) 

cs 
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Let  u(r,  <j>)  be  the  horizontal  velocity  in  the  boundary  layer  at  the  end  of  the  cylinder. 
Because  the  tow  direction  is  at  an  angle  a  to  the  cylinder  axis,  the  velocity  at  the  inlet  to  the 
control  volume  is  U0  cos  a.  There  is  no  flow  across  control  surface  4  (see  figure  7),  so  evaluating 
equation  (18)  gives: 

(U0  -  u)Ax  cos  a  =  (v  •A\.  (20) 

Using  equation  (20),  equation  (19)  becomes 


p(U0  cos af 

where  dAx  =  rdrd<j) .  For  a  cylindrical  geometry,  the  momentum  thickness  is  defined  by  Cipolla 
and  Keith3  as 


Un 


1- 


u 


U, 


Oj 


\dA, 


(21) 


a+S 

G2  +  2ad  =  2  | 


U, 


o  \ 


u. 


oj 


Vdr. 


(22) 


Equating  the  right  sides  of  equations  (21)  and  (22)  and  then  simplifying  yields 


S 

xp{U0  cos  af 


=  G2 +2  aO. 


(23) 


Since  the  left  side  of  equation  (23)  consists  of  measured  quantities,  this  expression  enables 
evaluation  of  the  momentum  thickness  at  the  end  of  the  cylinder. 
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3.  RESULTS  OF  EXPERIMENTS 


3.1  WALL  SHEAR  FORCE  AND  MOMENTUM  THICKNESS 

The  total  measured  drag  Fx  as  a  function  of  tow  speed  U0  and  angle  of  attack  a  are  shown 
in  figure  8  for  the  four  different  materials  tested.  These  data  were  used  to  calculate  the  axial  wall 
shear  force  S  for  the  four  cylinders  using  equation  (12). 


Figure  8.  Total  Measured  Drag  Versus  Tow  Speed  and  Angle  of  Attack 


Figures  9a  and  9b  show  the  data  plotted  as  S  versus  a  for  constant  speeds  and  S  versus  U0 
for  constant  angles  of  attack.  As  expected,  the  shear  force  for  a  fixed  angle  increased  with  tow 
speed;  however,  no  clear  trend  was  observed  when  the  angle  was  varied  at  a  fixed  speed. 


Tow  Angle  (°) 

Figure  9a.  Shear  Force  as  a  Function  of  Tow  Angle 
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Figure  9b.  Shear  Force  as  a  Function  of  Tow  Speed 


The  axial  wall  shear  stress  rw  and  the  tangential  drag  coefficient  Cd  were  calculated  from 
the  shear  force  S  using  equations  (13)  and  (14).  In  defining  Cd ,  the  tow  speed  U0,  rather  than  the 
axial  component  U0  cos  a,  was  used.  Figure  10  shows  r„.  versus  a  for  all  combinations  of  tow 
speed  and  angle.  Recall  that  the  wall  shear  stress  has  been  averaged  over  the  entire  surface  area 
of  the  cylinder  (circumference  and  length),  so  that  details  of  the  boundary  layer  asymmetry  are 
not  identified.  The  different  colors  denote  angles  that  are  similar,  and  the  symbols  denote 
equivalent  tow  speeds. 
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7° 
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= 

8° 
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Figure  10.  Mean  Wall  Shear  Stress  as  a  Function  of  Tow  Angle 
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Nondimensionalizing  rw  to  obtain  Q,  as  shown  in  figure  11,  leads  to  a  clear  trend  of 
decreasing  values  of  Cd  with  increasing  speed  over  the  range  of  angles.  Each  data  set  is  for 
constant  speed,  except  for  the  data  obtained  with  the  6.10-m-long  cylinder,  which  was  towed  at 
three  different  speeds.  A  large  variation  in  Q  was  measured  depending  on  the  value  of  U0. 
Also,  at  a  given  speed,  Q  increases  with  increasing  angle.  The  values  of  momentum  thickness 
0{L)  at  the  end  of  the  cylinder  were  calculated  using  equation  (20)  and  are  shown  as  a  function 

of  tow  angle  in  figure  12. 


Figure  11.  Tangential  Drag  Coefficient  as  a  Function  of  Tow  Angle 


Figure  12.  Momentum  Thickness  at  the  End  of  the  Cylinder  as  a  Function  of  Tow  Angle 
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The  expected  trends  of  increasing  6(L)  with  decreasing  tow  speed  and  increasing  angle 
are  apparent.  This  trend  indicates  that,  for  a  given  speed,  increasing  the  tow  angle  leads  to  the 
development  of  a  thicker  turbulent  boundary  layer.  The  dependence  of  wall  shear  stress  on  angle 
for  a  fixed  speed  is  somewhat  more  complicated  (as  shown  in  figure  10).  The  trend  of 
decreasing  mean  wall  shear  stress  that  would  be  associated  with  a  thicker  flat-plate, 
zero-pressure-gradient  turbulent  boundary  layer  is  not  apparent  and  suggests  a  departure  from 
that  ideal  case.  The  momentum  thickness  values  provide  an  outer-boundary-layer  length  scale 
for  subsequent  data  analysis. 

To  obtain  a  collapse  of  the  data,  the  nondimensional  quantity  (LI Of  was  introduced  as  a 
scaling  on  C</,  where  n  is  an  empirically  determined  constant.  This  scaling  introduces  the  only 
measured  boundary  layer  outer-length  scale,  with  the  ratio  LI 6  serving  as  an  indication  of 
boundary  layer  growth  over  the  length  of  the  cylinder.  The  results  are  shown  in  figure  13,  with 
n  =  0.5.  Included  in  these  results  are  a  limited  number  of  measurements  obtained  for  a  6.1-m 
cylinder;  these  measurements  collapse  with  the  3.05-m  cylinder  data.  It  is  interesting  to  note  that 
the  scaling  leads  to  a  linear  trend,  approaching  a  nearly  constant  value  of  Q  (LtOf.  Because 
each  quantity  was  varied  independently,  this  collapse  demonstrates  that  LIO  is  an  appropriate 
scaling.  In  addition,  it  is  shown  that  the  cylinder  lengths  are  sufficient  to  achieve  statistically 
meaningful  boundary  layer  flows  in  terms  of  the  spatial  development. 


Figure  13.  Nondimensional  Scaling  for  Drag  Coefficient 
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The  data  are  replotted  in  figure  14  for  data  sets  of  approximately  constant  values  of  angle 
of  attack,  with  Cd  plotted  as  a  function  of  U0  for  all  angles  obtained  for  the  3.05-m-long 
cylinders.  In  this  form,  the  data  display  the  expected  trend  of  decreasing  Q  with  increasing 
speed  for  all  angles — although  the  values  are  lower  than  those  that  flat-plate  theory  would 
predict  at  the  highest  speeds.  The  data  for  different  angles  follow  a  similar  trend,  with  the  largest 
angles  corresponding  to  higher  values  of  Q.  Figure  1 5  is  a  plot  of  uT  /U0  as  a  function  of 

Rex,  with  data  from  Coles7  shown  as  a  solid  blue  line.  In  this  Reynolds  number  range, 
corresponding  to  2  x  104  <  Ree  <  5  x  104,  the  quantity  uT/U0  decreases  with  increasing  length 
Reynolds  number  for  both  the  cylinder  and  flat-plate  data.  The  cylinder  data  highlight  the 
variation  of  shear  stress  with  small  changes  in  angle  at  all  speeds. 


Speed  (m/s) 


Figure  14.  Tangential  Drag  Coefficient  as  a  Function  of  Tow  Speed 


Figure  15.  Nondimensional  Friction  Velocity  as  a  Function  of  Reynolds  Number 
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Although  the  boundary  layers  are  asymmetric  for  cases  of  nonzero  angle,  they  can  be 
characterized  by  circumferentially  averaged  inner-  and  outer-length  scales.  The  ratio  of 
outer-length  to  inner-length  scale  values  (&v/uT)  is  plotted  as  a  function  of  Reg  in  figure  16, 
along  with  values  for  a  flat-plate  boundary  layer  taken  from  Coles.7  The  measured  values  do  not 
deviate  significantly  from  those  for  a  flat  plate,  and  no  clear  trend  is  apparent  with  speed  or 
angle.  This  result  is  expected  in  terms  of  the  ratio  S/d.  Assuming  that  8  is  an  order  of 
magnitude  greater  than  0,  the  ratio  S/d  is  then  of  order  1,  which  is  not  sufficient  for  strong 
curvature  effects.  In  contrast,  measurements  obtained  for  higher  aspect  ratio  cylinders3'5  showed 
values  of  6/  v/ur  that  were  two  to  three  times  higher  than  flat-plate  values  in  this  range  of 
Reynolds  numbers. 
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Figure  16.  Ratio  of  Outer-  to  Inner-Length  Scales  as  a  Function  of  Reynolds  Number 


Figure  17  shows  rw  versus  length  Reynolds  number  ReL  for  all  combinations  of  tow  speed 
and  angle  for  the  3.05-m  cylinders.  Also  plotted  as  open  circles  are  the  values  of  the  spatially 
averaged  r„  predicted  for  a  flat-plate  turbulent  boundary  layer8  of  the  same  length  at  zero  angle 
of  attack.  Consideration  of  the  variation  of  r„.  as  a  function  of  UoLJv  facilitates  comparison  of 
the  measured  data  with  flat-plate  boundary  layer  results.  The  length  scale  L0  can  be  viewed  as  a 
reference  length  (since  data  for  only  one  length  of  cylinder  are  included)  and  increasing  values  of 
UoLJvoxt  strictly  the  result  of  an  increase  in  tow  speed.  At  low  speeds,  the  values  of  r„  agree 
with  the  flat-plate  values;  however,  as  U0  is  increased,  r„  appears  to  reach  a  maximum  value  and 
departs  significantly  from  the  flat-plate  values.  This  result  corresponds  to  the  lower  Cd  values 
shown  in  figure  12  at  high  speeds.  This  behavior  was  observed  for  all  tested  angles  of  attack, 
including  the  smallest  value  a=  1.6°. 
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Figure  1 7.  Mean  Wall  Shear  Stress  as  a  Function  of  Length  Reynolds  Number 


If  it  is  assumed  that  a  well-defined  turbulent  boundary  layer  exists,  qualitative  aspects  of 
the  hairpin  vortices  in  the  boundary  layer  can  be  considered.  For  a  flat-plate  boundary  layer,  as 
U0  increases,  the  wall  shear  stress  increases  and  the  hairpin  diameter  and  spacing  decreases.  For 
the  cases  investigated  here,  the  hairpin  diameter  and  spacing  reached  a  minimum  and  then 
increased  as  U0  increased.  Assuming  an  average  hairpin  diameter  of  10  viscous  lengths  and  an 
average  spacing  of  100  viscous  lengths,9  the  maximum  number  of  hairpin  pairs  around  the 
circumference  of  the  cylinder  is  200,  based  on  the  maximum  value  of  zw.  Because  the  hairpin 
spacing  is  two  orders  of  magnitude  smaller  than  the  diameter  of  the  cylinder,  the  effects  of 
curvature  are  expected  to  be  small.  Therefore,  the  dominant  effect  is  likely  to  be  the  crossflow 
component  caused  by  the  nonzero  angle  of  attack.  Additional  experiments  in  which  cylinder 
diameter  is  varied  would  clarify  whether  the  value  of  200  hairpin  pairs  is  in  fact  a  meaningful 
constant.  Detailed  turbulence  measurements  would  reveal  if  the  flat-plate  assumption  of  hairpin 
spacing  is  valid,  or  whether  the  entire  structure  changes  for  this  class  of  boundary  layers. 


3.2  ORTHOGONAL  FORCE 

The  orthogonal  drag  coefficient  was  calculated  from  the  orthogonal  force  using 
equation  (17).  Figure  18,  which  is  a  plot  of  C„  versus  Rep  for  the  3.05-  and  6.1-m  cylinders, 
shows  that  there  is  a  significant  deviation  from  the  classical  cylinder  in  the  crossflow  case.  In 
this  case,  the  velocity  scale  is  taken  as  U0  sina,  which  is  appropriate  for  crossflow.  Each  curve 
represents  data  for  a  given  cylinder  density.  In  this  range  of  Reynolds  numbers  based  on 
diameter  (for  example,  4000  <  Rep  <  20,000),  the  drag  coefficient  value  for  pure  crossflow 
ranges  from  1.0  to  1.2  (as  published  by  White10),  which  is  approximately  2  to  10  times  the  values 
measured  in  this  experiment.  Not  only  are  the  values  significantly  different,  but  there  is  a  strong 
dependence  on  the  tow  angle,  which  is  inherent  in  the  values  of  C„  and  Rep,  by  virtue  of  the 
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sin(a)  term  in  each  quantity.  Figure  19  shows  the  orthogonal  drag  Fr  as  a  function  of  tow  speed 
for  approximately  constant  angles.  At  a  given  speed,  the  orthogonal  force  on  the  3.05-m 
cylinders  increases  with  increasing  angle,  as  expected;  however,  the  values  of  Fr  for  the  6.1-m 
cylinder  do  not  vary  over  the  range  of  angles  4.8  to  6. 1  °. 


Figure  18.  Orthogonal  Drag  Coefficient  as  a  Function  of  Diameter  Reynolds  Number 


Figure  19.  Orthogonal  Drag  Force  as  a  Function  of  Tow  Speed 
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To  collapse  the  data,  scaling  is  sought  for  the  orthogonal  force  that  takes  into  account  the 
effect  of  the  turbulent  boundary  layer  in  the  axial  direction.  From  the  data  in  figure  18,  it  is 
apparent  that  the  effect  of  the  boundary  layer  is  to  reduce  the  momentum  deficit  on  the 
low-pressure  side  (upper  surface)  of  the  cylinders,  resulting  in  a  lower  crossflow  drag.  This 
effect  is  likely  accompanied  by  a  marked  change  in  the  vortex  shedding  with  respect  to  the  case 
of  simple  flow  past  a  cylinder  without  the  presence  of  an  axial  boundary  layer.  An  inner  variable 
scaling  for  the  boundary  layer  is  taken  to  reflect  the  changes  in  the  crossflow  vortex  shedding; 
therefore,  the  force  Fy  is  normalized  by  the  corresponding  measured  axial  shear  force  values,  the 
results  of  which  are  shown  in  figure  20.  The  term  U0  sin«/wr  represents  a  Reynolds  number 
based  on  the  viscous  length  v/ur.  The  results  show  that  this  scaling  collapses  all  data  over  the 
range  of  Reynolds  numbers,  tow  angles,  and  lengths.  Although  additional  measurements  are 
required  to  refine  this  scaling,  this  analysis  provides  a  simple  method  to  estimate  the  orthogonal 
force  from  the  tension  (average  shear  force)  and  local  angle.  This  approach  can  be  used  to 
improve  models  of  array  shape  that  currently  use  a  single  average  value  for  the  crossflow  drag 
coefficient  for  all  angles  and  speeds. 


Figure  20.  Orthogonal  Drag  Force  Normalized  with  Shear  Force 
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4.  SUMMARY 


Direct  drag  measurements  on  long,  small-diameter  towed  cylinders  at  small  angles  of 
attack  and  high  Reynolds  numbers  were  performed  in  the  High-Speed  Towing  Basin  at  the  Naval 
Surface  Warfare  Center.  The  primary  purpose  of  this  investigation  was  to  determine  the 
sensitivity  of  the  tangential  drag  force  to  small  angles  for  a  long,  slender  cylinder  under  tow. 

The  use  of  a  high-speed  tow  tank  enables  development  of  very  high  Reynolds  number  turbulent 
boundary  layers  without  interference  from  the  walls.  The  inherent  design  of  the  experiment 
precluded  an  independent  variation  of  angle  and  tow  speed.  The  facility  itself  limited  the 
maximum  tow  angle  and  the  minimum  Reynolds  numbers  that  could  be  achieved.  Nevertheless, 
a  wide  range  of  parameter  values  was  covered  during  these  experiments.  The  temporal  and 
spatial  mean  wall  shear  stress  along  the  cylinder  surface  was  computed  from  the  measured  drag, 
and  the  data  deviated  significantly  from  flat-plate  boundary  layer  data.  The  tangential  drag 
coefficient  was  found  to  be  very  sensitive  to  small  changes  in  tow  angle  at  a  given  speed.  The 
ratio  of  cylinder  length  to  maximum  momentum  thickness  was  used  to  collapse  the  data  for  all 
angles  and  tow  speeds. 

As  tow  speed  was  increased,  the  value  of  tw  increased  to  a  maximum  and  then  decreased 
at  higher  speeds  for  all  angles.  This  trend  points  to  a  true  Reynolds  number  effect,  and  detailed 
turbulence  measurements  are  required  to  determine  the  alteration  of  the  boundary  layer  structure. 
The  trend  of  increasing  momentum  thickness  6{L )  with  decreasing  tow  speed  and  increasing 
angle  was  found.  At  a  fixed  speed,  increasing  the  tow  angle  led  to  the  development  of  a  thicker 
turbulent  boundary  layer.  From  these  measurements,  the  detailed  structure  of  the  asymmetric 
boundary  layer  are  unknown,  as  are  the  variation  of  rw  and  #  along  the  length  of  the  cylinder. 

The  force  orthogonal  to  the  cylinder  was  also  computed  from  the  data  and  was  found  to  be 
sensitive  to  tow  angle  and  speed.  The  related  orthogonal  drag  coefficient  is  significantly  lower 
than  values  for  cylinders  in  pure  crossflow.  The  measured  axial  shear  force  was  used  to 
normalize  the  orthogonal  force  and  to  collapse  the  data  for  all  angles  and  tow  speeds.  Additional 
experiments  using  cylinders  of  additional  lengths  and  diameters  are  required  to  further  refine  the 
scaling  laws  for  orthogonal  force  and  tangential  drag  coefficient. 

The  results  of  these  experiments  confirm  that  simple  corrections  to  drag  forces  using 
components  of  flow  based  on  geometry  are  insufficient.  The  sensitivity  of  cylindrical  boundary 
layers  to  small  changes  in  angle  of  attack  agrees  with  the  findings  of  Willmarth  et  al.1  For  towed 
array  shape  modeling,  these  results  show  that  applying  flat-plate  theory  to  determine  the  values 
of  the  tangential  drag  coefficients  will  lead  to  significant  errors.  Assumptions  for  the  orthogonal 
drag  based  on  simple  crossflow,  without  the  presence  of  a  boundary  layer,  are  also  inadequate. 
While  the  geometries  investigated  in  this  research  represent  a  small  subset  of  those  encountered 
in  towed  array  maneuvering,  they  establish  a  set  of  baseline  measurements  for  this  class  of 
problems.  Future  studies  should  address  additional  values  of  tow  angle  and  non-zero  yaw 
angles,  with  the  focus  on  the  effect  of  the  turbulent  boundary  layer  on  this  class  of  flows. 
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